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ABSTRACT 


This  report  describes  the  continued  evaluation  of  the 
Norwegian  short-period  Seismic  Array  (NORSAR),  which  was  conducted 
by  Texas  Instruments  Incorporated  at  the  Seismic  Data  Analysis  Center 
over  the  period  l  April  1972  to  31  March  1973. 

The  major  areas  of  study  presented  in  this  report  are: 

•  Signal  analysis  on  a  regional  basis 

•  Signal  amplitude  response  patterns  across  the  array 

•  Array  processing  performance 

•  NORSAR  seismic  event  detection  capability 

•  Behavior  of  short-period  seismic  discriminants 

The  total  data  base  for  this  study  comprises  344  events. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
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SECTION  I 
INTRODUCTION 


This  report  presents  the  results  of  an  ongoing  evaluation  of 
the  short-period  (SP)  Norwegian  Seismic  Array  (NORSAR),  using  seismic 
data  recorded  during  1971  and  1972.  The  overall  objectives  of  the  NORSAR 
SP  evaluation  are: 

•  Determine  the  best  processing  methods  for  enchancing  the 
signal-to -noise  ratio  of  Eurasian  events. 

•  Determine  the  array  capability  for  Eurasian  events  . 

•  Evaluate  the  performance  of  short-period  discriminants 
at  NORSAR. 

•  In  conjunction  with  long-period  NORSAR  data,  determine  the 
detection  and  discrimination  capability  of  NORSAR  for  Eurasian 
events . 

Substantial  progress  has  been  made  toward  achieving  the  first 
three  objectives,  but  the  results  presented  in  this  report  may  still  be  improved 
as  the  data  base  for  the  evaluation  is  expanded.  Work  toward  meeting  the  fourth 
objective  is  still  at  an  initial  stage,  and  will  be  given  more  attention  in  our  fu¬ 
ture  analysis. 

Five  analysis  tasks  were  undertaken  in  order  to  meet  the  first 
three  objectives  stated  above: 

0  Noi .  j  analysis 
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Signal  analysis 

Array  processing  effectiveness 
Detection  threshold  estimation 
Behavior  of  SP  discriminants 


Results  from  noise  analysis  were  presented  in  Special  Report 
No.6  under  the  Extended  Array  Evaluation  program.  No  further  vork  on  this 
subject  has  been  undertaken  since  then.  Results  from  each  of  the  four  remain 
ing  tasks  are  presented  in  subsequent  sections  of  this  report. 


The  NORSAR  SP  array,  centered  about  100  km  due  north  of 
Oslo,  Norway,  consists  of  132  short-period  seismometers  and  has  an  aper¬ 
ture  of  about  100  km.  The  sensors  are  grouped  in  22  six-element  subarrays; 
each  subarray  has  a  center  sensor  and  a  five-sensor  ring  and  is  about  7  km 
in  diameter  (Figure  1-1).  Throughout  this  report,  the  official  nomenclature 
will  be  used  whenever  a  NORSAR  seismometer  or  subarray  is  referred  to. 


The  results  presented  in  the  following  sections  are  based  pri¬ 
marily  on  events  located  on  the  Eurasian  Continent.  A  total  of  344  events 
have  been  analyzed;  106  of  these  were  listed  as  the  data  base  for  Special  Re¬ 
port  No.  6.  Table  1-1  lists  the  parameters  for  the  additional  238  events  pro¬ 
cessed  since  then. 


Geographically,  the  events  are  concentrated  principally  in  the 
Northwestern  Pacific  (from  Kamchatka  to  Taiwan)  and  in  south  central  Asia 
(north  and  west  of  the  Himalayan  system).  Thirty-nine  events  from  the 
Mediterranean  region  are  included  as  well  eight  from  the  Arctic  Ocean  and 
eight  from  Continental  North  America.  Twenty-one  events  are  presumed 
explosions;  including  eight  from  Eastern  Kazakh,  one  from  Western  Kazakh, 
two  from  the  Ural  Mountains,  three  from  Western  Russia,  one  from  Novaya 
Zemlya,  one  from  the  Aleutian  Islands  and  five  from  Nevada. 
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70. 6T 

NF  P 

4.0 

R 

04/14/71 

14.75.94 

96. 7  N 

l11  .5F 

vnt> 

4.4 

D 

04/19/71 

l  4. 04. C 9 

9  ?  .  0  M 

140. sr 

9  9 

6.1 

r> 

04/19/91 

77.17.39 

4  l  .  6  \ 

79. 2  p 

MOP 

4.9 

D 

04/17/71 

09.77.  ri9 

4  4 . 4  u 

1  4  3 . 9  P 

NIP 

4.9 

R 

04/10/71 

I  7.7  3.02 

41 . 9  \ 

70.9c 

NOP 

r,  n 
'  •  / 

0 

04/10/71 

21. 09. 49 

4  1  .  9  M 

70.4c 

N'P 

4.7 

r> 

07/00/71 

14.44.19 

43. 9  \ 

147. 79 

4  6 

4.9 

r> 

07/10/71 

03.09.  00 

4  3.6  N 

147.79 

36 

4 . 3 

p 

07/10/71 

09. 01. 74 

4  9  .  ON 

19  0.99 

Nn« 

4.6 

p 

07/12/71 

02.12. 70 

93.1  M 

160.0P 

VPP 

4.o 

o 

09/16/71. 

04. 9P. 00 

7  p  .  0  N 

103. 7F 

MOP 

9.9 

P 

o p/14/71 

1 3.7C. 24 

2  9 . 3  A1 

1 03. 7 F 

NOP 

'•.9 

P 

0  9/19/71 

1 4. CO, 90 

37.  1  N 

1 1  6  .  ^  W 

0 

9.4 

P  F 

0  °  /  l  °  /  7 1 

09.29. 97 

2  7.  QV 

1  7  1 . 9  c 

-i  0 

t  . 

9.4 

P 

09/1 °/71 

22.  lc. 37 

4  9 . 3  M 

199. 4F 

VHP 

6.0 

D 

0H/?7/71 

1  7. 94.  14 

30.1  N 

60.79 

NIP 

9.1 

P 

0°/?.?/71 

09.36.  11 

7  7 . 0  A' 

107.79 

NOP 

9.7 

D 

09/27/71 

21.99.17 

4  9  .  (■  N 

1  9  1 . 0  c 

74 

9.7 

P 

09/24/71 

16.33. 77 

5  2.2  M 

91  .49 

MOR 

5.7 

0 

09/79/71 

00.30.44 

2  '3  .  2  N 

9  7.30 

Nn» 

4.1 

P 

09/29/71 

10.41.4? 

4  7  .  1 N 

147.79 

6  9 

4.1 

P 

09  / 24/71 

06. 95. OP 

70. ON 

c  0. 7C 

49 

4.3 

P 

09/27/71 

09.2C.19 

30.  7N 

9Q.  7P 

9  4 

9.0 

P 

09/27/71 

07.99.11 

30.  IN 

90.  7P 

63 

4.4 

P 

09/77/71 

1 3.49.59 

4  0.  «  N 

1  4  7 . 3  C 

NOP 

4.9 

P 

09/29/71 

1  9. 9  7.  4  7 

2  9 . 3  N 

1  30.79 

39 

9.7 

P 

09/79/71 

16.74.44 

3  7.6N 

99.  wc 

NOP 

4.3 

P 

Q9 /2C /7l 

15.16.56 

7  4.9  K 

7  3 . 9  f- 

NfP 

9.0 

D 

0  9  /  •*  i  /  7 1 

01. 92.  1  6 

36. 9N 

7M.RO 

,3  7 

4.7 

n 

00/02/71 

09.70. 27 

3  .  n  N 

r  *.7r 

345 

4.7 

P 

00/07/71 

19.74.47 

3  0.  1  M 

90.3c 

49 

5.1 

P 

00/02/71 

77.7i.39 

7  0.  1  N 

9  0 . 3  F 

79 

5.0 

P 

00/03/71 

19.42.16 

2  3.9V 

103.7F 

NOP 

4.3 

D 

09/09/71 

12. 1 9.96 

37. 7AJ 

30.  3T 

20 

4.9 

P 

09/09/71 

21.1 3.57 

4  b  .  6  N 

1  4  0 . 9  C 

N'T  D 

4.6 

P 

09/04/71 

06.49. 99 

46.4  V 

14  1.1  C 

16 

5.7 

P 

09/04/71 

13.37.  10 

46. 7N 

1  4  l  .  4  F 

7  c 

6.1 

P 

09/04/71 

20.  10.4.7 

4  6 . 6  N 

1 4  1 . 2  C 

MOP 

5.0 

P 

09/09/71 

04.10.  19 

49. 9N 

27. OF 

140 

3.4 

P 

00/09/71 

O'7. 29.  14 

?7. ?N 

1 4  1 . 7  c 

96 

9.6 

P 

09/09/71 

11 .4P.27 

4  6.4V 

14 ) .4r 

6 

5.9 

n 

09/09/71 

1 7. 97.34 

2  9 . 7  m 

6  O.OP 

1  7 

5.4 

P 
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rv  r 

nr s  ion  at  j  nv 

FPS/251/16M 
TUO/2S1/17M 
n  S  /  7  ^  1  /  1  o^| 
T1  IP  /  ?  5  1  /  ??v 
p  R  S  /  ?  5  ?  /  ?  1 M 
KI.IR/25?/73\) 
A^n/p  ms/oSM 
HIV/  ’  56/Q1M 
K  l  JP  /?5A/  1  ?»| 
Hr'\l  /  P  5  H  /  1  4N 
R  AT/? A  1  /07V 
VFV/  ?6?/OON 
WFW96?/  1  1M 
2 64/0  IN 
T  F°/?6A/ot>M 
T'JP/764/1 AM 
FR  S/7,65/  ]  4\' 

AFG/?68/0«m 
0  8  F  /  ?  a  r,  /  o  5  N 
TAL/269/  15V 
HIV/26?/? 3M 
TIJR/271/05V 
TAL/271/1  AN 
MnM/271/1AM 
MFV/272/14M 
TIIR  /  2  7  3  /  n  PN 

p  AT/273/  l  IN 
ST  P/27  3/ ? jm 
TIJR  /?  7f  / 0 7N 
TIJO/276/17N 
SI R/278/niM 
IP  A/27R/1RN 
TMR/27R/IRN 
TR  A/27R/72N 
TIIR/27Q/01N 

PAT/279/  1  IN 
NPV/29  1/1  a\i 
KA7/282/06N 
I  R  A/?°R/ 1 4M 
AFO/288/ L6M 
K  A  7  /  3  5  6  /  0  6  N 
S  ID/0Q1/15N 
KIIP/001/16V 
K'lR/DOl  / 1  Rm 
K * .1 P  /002/05N 


n  A  T  R 

n^/oR/7] 
n  R  /  O  R  /  7 1 
n  o  /  o  0  /  7 1 

OP/qp/71 
OC  /OR/7] 
OG/09  n 1 
n  7  / 1  7  /  7 1 
OF/ I  3/71 
OR/ 1 7/ 71 
07/15/71 
00/18/7] 
07/17/71 
OF/10/71 
09/21/7] 
0  c  /  2  1  /  7 1 
09/21/7] 
00/00/7] 

0°/?  S/71 
07/25/7] 

0F/R5/71 

0  c  /  2  5  /  7  1 

09/28/71 

07/28/71 

OR/28/71 

0n/?7/71 

09/20/71 

09/^0/71 

09/70/71 

10/03/71 

10/07/71 

10/05/71 

10/05/71 

10/05/71 

10/08/71 

10/06/71 

10/06/71 

lO/On/7] 

1 0/0r / 7 1 
10/15/7] 
10/15/71 
12/22/71 
01/0! / 7? 
01/01/7? 
01/01/7? 

01/07/7? 


no  jr  j  n 
T  T  Mr 

1  6  .  5  c.  .  a  p 

1  7.  01. 09 
19.7?. 15 

2  2.75.  15 
21.06.20 
2  7.01 . 06 
Oc. Q7. ?5 
01 , 45.  35 

1 a.01.2o 

I  4  .  5  5.  Q  C 
02.12.7c 
00. 58.  35 
1  1 . 00. OA 
01  .  C4.  1  7 

0  n .  n .  5 1 

16.  40.  51 

1  A.  20.  10 

0  3.  57.00 
05.44. 34 

1  c  .  2  F  .  1  p 

2  7.4  8.  2 5 
C5.  1 C. 2 A 
14.  04.41 
1  a.  n.  Oo 

1A. 00. 00 

08.55.50 
11.5?.  36 
2 1  .  3 1 . 2  5 
07.  aa.  ?a 
1  7.  ]  P.  57 
01  . AO. 41 
18.71.17 

1  8.  63. 04 
77. AS. 04 
01 . 46.  7  p 
11.77.  44 
1  A.  30.  00 
06.02. 57 
1  A.  ]  <*.  71 
l 6. 22.  17 
06. 5F. 56 
1 5.04.  19 
1  6. 55. 06 
1  8.  13.64 

05.7  7.  pc 


p  pij  pr  = 


46. 6* 

44. 4  N 
3  7.0V 
75.  7N 
A  R  .  5  N 
3  9.1  N 
51  .ON 
7  9.4V 
5  7  .  «  V 

7  8.3V 
3  2  .  A  M 
37. 3N 

46.4  V 
37.7  V 
37. ov 

7.2  V 
7  6.3  M 
37.  1  V 
3.  7 V' 
A  0.  3  A 

3  7.0V 
3  7.7V 
51.7V 
61  ,6  V 
3  8 .  o  V 
36.7V 
6’T.  s  v 

2  7.7V 

3  9  .  O  N 

3  1.6V 

3  8.7V 
52.1V 
77.  1  V 
5  0.0V 
7  7.3V 
37.9V 

4  7  .  a  V 

5  9.7V 
50.7V 
A  9.4  V 
A6.  IV 


1  nv 

RFPTH  w 3 

RLTV 

M  160.9- 

16 

5.9 

D 

J  7  0 . 3  r 

6 

A.  R 

P 

l  1  A  0 . 9  F 

1  9 

6.3 

P 

,1 

17 

4. 8 

P 

1  A  0  .  o  p 

MOP 

4.5 

P 

i  1  6  0  .  n  p 

7 

6.0 

P 

•  7 1 .  7  r 

1  10 

4.5 

P 

1  A  9 . 0  c 

1  ?  9 

4.0 

n 

1  14  8.91- 

36R 

4.6 

P 

1  1  A  7 . 4  F 

17 

5.0 

D 

1 7p. 6P 

1  1? 

A  .6 

P 

11 8. OF 

10 

A. 6 

P 

A]  .  1  r 

NOP 

A  .5 

P 

A-R.cn 

19 

6.2 

P 

cl  .op 

MOP 

6.0 

P 

3  0 . 3  F 

37 

A.  7 

P 

1  a  n .  3  p 

14 

A. 7 

P 

A  9 . 7  p 

56 

A.  5 

P 

27.  ?F 

107 

A  .  1 

D 

1  ?  5 . 7  c 

67 

5.3 

P 

70.29 

2  2  9 

A  •  6 

P 

30. 1  c 

36 

a. 7 

P 

126. OF 

169 

5.9 

P 

1  A  7 . 6  C 

45 

A.  2 

P 

1  1  6 . 0  !ri 

0 

A  .4 

D 

30.  P 

NOR 

A. 5 

P 

177. pc 

A  1 

5.0 

P 

1 A  n .  3  c 

NOR 

5.6 

O 

79. or 

7  3 

6.7 

P 

79.99 

27 

6.6 

P 

1  7  2 . 8  U 

vnp 

5.3 

P 

5  6  .  a  r 

7  Q 

5.1 

P 

79. 7F 

9 

4.5 

n 

C0.7L‘ 

6  6 

4.1 

P 

30.2  8 

19 

4.6 

p 

1.78. 2C 

1  39 

4.9 

D 

1 1 1> .  0  W 

0 

A  .  7 

D 

77. 7E 

0 

6.6 

P 

54 .69 

7  c 

A  .  7 

P 

71.0" 

167 

A. 9 

P 

48. ?F 

0 

6.0 

P 

1 5? .pr 

NOP 

4.1 

S 

156. DP 

MPO 

A. 6 

7 

1 56.59 

50 

4.0 

s 

1A6.?P 

NOR 

4.0 

s 

m  n  r 
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FVFNT 

DFS  IGNAT  J-)M 

PAT  r 

PR  1 0 1  N 

T  I  M- 

l  AT 

LON 

nrpi  H 

MO 

s  nu  p  r.  8 

01  TN 

OR8/002/09N 

0  1  /O  2  /  72 

0°. 17. 53 

3  7  .  c  n 

20. 7F 

NOP 

4.? 

s 

S  IM/00?/10»! 

01  /O  7 / 7  2 

10. ?  7. 35 

41 

0  4 . 5  F 

1  c 

8.7 

c. 

K  A  M  /OO?  /06N 

01  /O'5/  7? 

04. 36. 3  R 

5  1 . 4  N 

1  rO,4F 

NOR 

4.8 

P 

KAM/OO^/IPM 

01/07/72 

19.24.  43 

52.  ON 

1  5  0 . 0  f 

Npc 

4.8 

N 

K AM/004/0PN 

01/04/72 

02.?9.  1  P 

5  5 . 6  N 

161 .2F 

NPP 

4.3 

c 

K  AM /004 / 1  ON 

01/04/7? 

10.42.71 

5  5.6M 

1  4  3  .  P  F 

NCR 

4.4 

S 

TAT/004/12M 

01 /04/72 

1  7. 1 5.  17 

2  2 . 4  N 

17?. 78 

Mf]  (> 

4.H 

P 

K(JR/005/0?M 

01/05/7? 

02.14.10 

47. OM 

147.? p 

Nfl« 

4.4 

p 

AUS/005/04M 

01/05/72 

0^  .57.  41 

4  7 . 0  N 

!  4.7P 

1  1 

4  .0 

P 

tao/oo5/i:n 

0  1  /Or'/ 7? 

1  ?.0?.54 

37.  ^ 

77. 1  r 

NP» 

4.8 

S 

KflM/005/1  4M 

01 /05/7? 

14.26.4P 

54. 6.N 

140. 4F 

NO” 

4.0 

5 

K  AM / 00 5 / 1 4M 

01 /08/7? 

1  4.  C9.  50 

5  7 . 3  N 

1 60.5  r 

<>8 

7.Q 

5 

K IP /006/06M 

01/04/72 

P4. 70. 76 

4  0.  7  N 

7  ?  .  4  r 

NPR 

4.7 

D 

TA  I  / 00 6 / OoN 

0! /06/7? 

04. 7  p. 74 

23.7  K 

1  ?  3 . 4  r 

N'PR 

4.7 

n 

I f A/ 004 /OON 

01 / 0  4 / 7  2 

05.41.33 

3  0.  3  N 

80.  8C 

NPR 

8.? 

P 

SWP/n07/?0N 

01/07/7? 

?C.  7?.  32 

44.1  N 

48. 1  r 

NPR 

4.7 

5 

K0M/009/07N 

01 / Or / 7? 

07.77. 04 

5  4  •  4  N 

144.48 

MUR 

7.6 

S 

K  AM / 00°-  / 1  4N 

01/00/7? 

1 4. or. 

8  5.  7  N 

167.60 

NPR 

4.7 

c. 

K 1  IP  /QC9/14N 

01 /OP/72 

14.4  7.46 

48.  IN 

148.48 

NPR 

7.P 

s 

o  H  I  /  0  1  P  /  0  rN 

01/10/7? 

05. ?3. 5? 

?0.  cjv 

}  70.48 

NPR 

5.0 

s 

KOM/01 I /08N 

01/11/72 

08. 54. 74 

54. 7M 

16«.?8 

7  P 

3.9 

s 

rpr/ni?/i3N 

01/12/7? 

17.51.20 

3  5 . 0  N 

?7 . 5  r 

NpP 

4.0 

0 

K AM/012/20N 

01/12/7? 

20.20.15 

55.  6N 

1  43.98 

NOR 

4  .  P 

S 

SI Q/01 3/17N 

Ol/n/7? 

17. ?4. 07 

41  .ON 

147. 1 F 

NOR 

5.7 

<; 

ST«/014/0?N 

01/14/7? 

07. ?c.  ?0 

67. 5  N' 

171 ,8F 

Nnt 

7.0 

p 

IP A/014/P2M 

01/1^/72 

7?. 1 0. 04 

3  ?  .  0  \ 

44.05 

N'op 

8.1 

p 

KUO/01  5 /09N 

01/15/7? 

00. 5F. 3? 

49.  6  M 

1 85. or 

NPR 

3.0 

S 

F«S'015/nN 

01/15/72 

!  p.  07.  5° 

5  7.4N 

120.78 

13 

4.7 

R 

S I  1/015/20N 

01/15/7? 

20. ?1. 50 

40. 3N 

79.08 

NPP 

5.4 

S 

M  2/01 5/20N 

01/15/7? 

20.45.  ?? 

3  9 . 7  N 

70.  or 

Nr;p 

4.6 

D 

K AM/016/04N 

01/14/72 

04 . 7o.  16 

55. 4N 

147.58 

NOR 

3.« 

s 

KAM/016/1 IN 

01/14/77 

1 1.00.4C 

5  5  .  ^  N 

1  4  7  .  r>  r 

?  5 

7.9 

s 

n A/01R/21N 

01 /l P/7? 

21.  12.0? 

3  7 . 5  \ 

4  R  .  7  F 

NpR 

A  .0 

p 

l TA/01R/27N 

01  /I  p/7? 

23. 14.12 

4  4 . 7  N 

p .  ?r 

->  R 

C. 

4.1 

p 

onn/o?o/02N 

01  / 20 / 7  ? 

07.  1  c. .  07 

3  4.4  N 

77. ir 

NOR 

4.P 

s 

KIJR/022/01N 

01/22/7? 

01.41. 24 

5  0.  ON 

I 8?. OP 

NPP 

4.7 

p 

TUP/022/1 7N 

01/22/7? 

1  7.  I  7.  7] 

3  7.  6,N 

?a.Q8 

17 

4.4 

p 

K  AM/Q?5/1 ON 

01/25/72 

10.02.40 

53.  ON 

1  4  0.  qc 

NPR 

4.6 

p 

IT A/025/20N 

01 /?4/7? 

20. 24.39 

4  3  .  P  N! 

1  7.48 

38 

4.5 

p 

fP  F/026/12N 

01/26/7? 

1?. 54. 39 

7  4.  5  N 

?  8 . 5  r 

N'PR 

4.0 

<; 

K AM/027/20N 

01/27/7? 

20.37.  ?P 

55. 7N 

142.38 

40 

3.8 

S 

FC S/02R/04N 

01/2P/7? 

04.??. 28 

2  7.5N 

126.88 

NPP 

4.4 

PAK/028/10N 

01/28/7? 

10.24.  54 

26. 6N 

66.3C 

NOR 

4  .0 

N 

KIP/02P/20N 

01  /?R/7? 

20.  2c.  1  9 

4  l.ON 

78.08 

NOR 

4.4 

P 

Fp S/OPR / ? IN 

01  /2  R / 7? 

21.50.00 

45. ON 

1 36. OE 

NOP 

4.0 

N 

rOW'H-NT 

«  pM 

MO 

M  ON] 

NO 

\'P 

MHN 

M  F>M 
mpi\j 
MP 

ruj 

M  R  N 

Mprj 


ttBMMB#ll*y#i  WBfcntSWUBl 

I 


»* 

fv^n  r 

pec;  IGNAT  I  ON 


Q 

0 

0 


... 

«  a 

«  *< 

Min 

I 


KIIR/028/23N 
I R  A/ 029/09N 
KAM/032/10N 
KAM/037/04N 
Kt)R/033/09N 
KUR/033/17M 
YtJN  /Q34/07N 
T  T A /03  6 /0?N 
BA  I /035/03N 
IT  A / 0  3  5 /C  4N 
ITA/O3^/09M 
TTA/035/17N 
TTA/035/19M 
ITA/036/niM 
IT A/036/03N 
TTA/036/05M 
ITA/03B/07N 
AOR/037/01N 
KA7/037/0qN 
ITA/037/21N 
I TA/039/12M 
IP  A / 04 1  /  R9N 
IR A/041/16M 
SIN/042/06N 
TIP/042/1  2M 
K  AM/Q42/21N 
K<)R  /044/05M 


I 

I 

I 

I 

I 

T 


r,R  p/044  /  I  3M 
K DM/044/ ??M 
KUP/046/16N 
OP  P/047/00M 
STN/047/23M 
KUP/049/1 PM 
F  IN/05  1/1 QN 
n  K  H  /  0  P  !  /  1  ON 
K  AM/Q6  1/20M 
K  AM/052/22N 
YU0/052/27M 
MOM/ 05 3/0 IN 
H I N /0  5  3/OPN 
KIJR/054/C3N 
K  A  M/054 /  19M 
KIJR/055/10M 
KUR/055/1PN 
KUR/056/19N 


*  0 
i 
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OA  T  P 
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S  :  SAAC/LASA  Bulletin 

P  :  PDE  Bulletin 

N  :  NTNF/NORSAR  Bulletin 

I  :  Bulletin  from  the  International  Seismic  Month  (ISM) 

ND  :  No  Detection  at  NORSAR  (TI  Analyst  decision). 

MBN  :  NORSAR  Magnitude 

E  :  Presumed  explosion 
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Data  quality  was  excellent;  about  one-half  of  the  time  all  132 
sensors  were  operational.  In  most  other  cases  one  subarray  (six  sensors) 
was  dead  or  contained  calibration  signals,  the  worst  data  loss  was  24  sensors. 
For  a  total  of  six  events  we  found  spikes  in  the  data,  but  these  events  could 
still  be  processed.  As  was  stated  in  Special  Report  No.  6,  it  appears  that  the 
seismometers  are  reasonably  well  equalized  across  the  array. 

From  about  March  20,  1971,  to  near  the  end  of  the  year,  at 
least  two  sensors  were  observed  to  have  reversed  polarities.  Sensor  06B03 
was  in  the  reversed  state  throughout  this  period,  and  sensor  05C04  was 
corrected  in  mid-August  only  to  be  followed  by  a  phase  reversal  of  08C02. 

For  about  ten  days  around  this  transition,  as  many  as  seven  sensors  were 
affected  by  this  problem.  No  phase  reversal  has  been  observed  in  1972  data. 
Table  1-2  shows  the  pattern  of  the  reversals  and  the  extent  of  the  data  base 
from  which  the  above  conclusions  were  determined.  The  effect  of  the  reversals 
on  subarray  and  array-beam  quality  will  be  discussed  in  Section  III. 
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SECTION  II 


SIGNAL  ANALYSIS 


A.  INTRODUCTION 

Several  large  signals  recorded  by  the  entire  array  were  analyz¬ 
ed  in  order  to  study  signal  characteristics.  Analyses  of  the  following  phenomena 
were  performed: 

•  Variation  as  a  function  of  epicentral  location  of  signals  in  both 
time  and  frequency  domains. 

•  Variation  of  SNR  (signal -to -noise  ratio)  from  subarray  to  sub¬ 
array  for  a  given  event,  and  its  relation  to  subsurface  velocity 
structure. 

•  Deviation  of  propagation  across  the  array  from  a  plane-wave 
model  and  the  variation  of  this  deviation  with  epicentral  loca¬ 
tion  and  oignal  spectral  content. 

•  Discrepancies  between  body-wave  magnitudes  (m^)  based  on 
measurements  of  data  processed  in  TI's  Alexandria  office  and 
those  reported  in  the  NORSAR,  LASA,  and  PDE  bulletins  and 
their  variation  with  epicentral  region  and  m^. 

Results  of  these  analyses  are  discussed  in  the  remainder  of  this  section. 

B.  REGIONAL  SIGNAL  CHARACTERISTICS 

A  study  was  undertaken  to  examine  signal  waveforms  recorded 
at  NORSAR  on  a  regional  basis.  For  this  purpose  it  was  found  convenient  to 


apply  the  short-period  discriminants  defined  in  Chapter  V  to  a  set  of  presumed 
earthquakes  from  Eurasia.  A  total  of  97  events  were  selected,  all  of  body  wave 
magnitudes  between  4.  5  and  5.0.  The  reason  for  selecting  this  narrow  range 
of  magnitudes  was  to  minimize  side  effects  on  signal  complexity  and  frequency 
contents  caused  by  magnitude  variations.  The  geographical  distribution  of 
these  97  events  is  shown  in  Figure  II- 1. 

F(  t  each  discriminant,  the  events  in  our  data  base  were  ranked 
according  to  their  discriminant  values.  The  25  highest  ranking  and  the  25 
lowest  ranking  events  were  then  plotted  in  geographical  space.  Figure  II-2 
shows  the  results  for  the  autocorrelation  mean  square  discriminant,  which  was 
selected  as  a  representative  for  the  discriminants  based  on  event  complexity. 
Similarly,  the  dominant  period  discriiv  inant  was  selected  among  the  criteria 
based  on  signal  frequency  contents,  an  1  the  resulting  plot  is  shown  in  Figure  II-3. 
It  should  be  noted  that  all  complexity  discriminants  showed  similar  results,  as 
was  also  the  case  with  the  frequency  discriminants. 

The  following  observations  on  regional  signal  characteristics 
may  be  made  from  Figures  II - 1 ,  II-2  and  II-3. 

•  Mediterranean  Region 

Not  unexpectedly,  signal  traces  from  this  close-in  region 
(around  20  degrees  distance  from  NORSAR)  tend  to  have  high 
complexity  values.  It  is,  however,  surprising  to  find  that 
most  of  our  events  from  the  Mediterranean  have  a  low  dominant 
frequency.  In  fact,  three  events  from  Turkey  and  one  from 
Italy  had  a  dominant  period  greater  tha*  1.  5  seconds  on  the 
adjusted-delay  array  beam,  and  all  18  events  from  Turkey/ 

Italy  exceeded  0.9  seconds. 
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Figure  II-4  shows  the  "average"  power  spectrum  for  our  Italy 
events  (i.  e  the  average  dB  value  for  the  selected  events  for 
each  frequency).  The  predominantly  low  frequency  content  is 
clearly  seen,  and  individual  inspection  of  the  spectra  shows 
furthermore  that  spectral  peaks  tend  to  repeat  from  one  event 
to  another.  This  is  not  so  surprising  in  view  of  the  proximity 
in  location  of  the  Italy  events. 

In  contrast  to  what  was  observed  for  Turkey  and  Italy,  events 
from  Greece  tend  to  show  significant  high  frequency  content. 

Four  of  the  five  Greece  events  included  in  our  data  base  had  a 
dominant  period  between  0.6  and  0.8  seconds  on  the  adjusted- 
delay  array  beam. 

Iran  /  Middle  East 

From  our  limited  data  base  it  appears  that  events  from  this 
region  generally  have  low  to  intermediate  complexity.  Dominant 
period  for  Iran  events  typically  ranges  from  0.6  to  0.9  seconds. 

Central  Asia 

Events  from  Central  Asia  have  in  general  low  complexity 
values.  This  is  especially  pronounced  for  the  Afghanistan/ 
Tadzhik/Sinkiang  region,  where  6  deep  events  were  included 
in  our  data  base. 

With  respect  to  dominant  period,  an  interesting  transition 
occurs  at  about  76  degrees  East  longitude.  Signals  from 
Afghani stan/ Tadzhik  west  of  this  meridian  are  generally  of 
high  frequencies,  with  typically  0.  5-0.8  seconds  dominant 
period  both  for  deep  and  shallow  earthquakes.  To  the  contrary, 


signals  from  the  Kir giz-Sinkiang  border  region  further  east 
exhibit  a  dominant  period  of  1.0-1.  3  seconds  on  the  adjusted- 
delay  array  beam. 

Signals  from  earthquakes  occurring  elsewhere  in  Central  Asia 
quite  often  indicate  a  significant  content  of  high  frequency  energy, 
but  our  data  base  is  still  to  scarce  to  permit  more  specific  con¬ 
clusions  for  this  general  area. 

•  Southeast  Asia 

Very  few  events  from  this  region  were  included  in  the  data  base, 
but  is  interesting  to  notice  the  high  complexity  of  our  Taiwan  ' 
events.  Signals  from  this  area  tend  to  have  a  dominant  low 
frequency. 

•  East  Asia  and  Western  Pacific 

Most  events  from  this  region  produce  low  complexity  signals  at 
NORSAR.  Some  high  complexity  signals  are  seen  for  events 
from  Kamchatka. 

There  is  a  considerable  variation  in  spectral  characteristics 
for  signals  from  the  Kamchatka -Kurile  arc.  The  two  extremes 
in  our  event  population  are  KUR/001/16N  and  KAM/078/18N 
with  dominant  periods  of  0.4  and  1.3  seconds  respectively 
for  the  adjusted-delay  array  beams.  Both  of  these  events 
are  of  unknown  depth.  In  general  events  from  this  region 
tend  to  have  a  dominant  high  frequency;  this  is  most  pronounced 
for  earthquakes  from  the  Kurile  Islands. 

In  conclusion,  there  seems  to  be  significant  variations  in  re¬ 
gional  signal  characteristics  observed  at  NORSAR.  Some  of  our  observations 
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contrast  with  the  general  expectations  of  decreasing  signal  complexity  and  in¬ 
creasing  dominant  period  as  the  epicentral  distance  increases. 

It  is  interesting  to  compare  the  regional  characteristics  found  in  the 
pi'eceeding  study  to  the  corresponding  characteristics  of  underground  explosions, 
which  are  low  complexity  and  low  dominant  period  (high  frequency).  The  most 
"explosionlike"  earthquakes  in  this  sense  are  seen  to  occur  in  the  Tadzhik/ 
Afghanistan  region  and  the  Kurile  Islands. 

C.  SUBARRAY  BEAM  AMPLITUDE  VARIATIONS 

The  striking  variations  in  signal  amplitudes  from  subarray  to  sub¬ 
array  were  discussed  briefly  in  Special  Report  No.  6.  A  study  was  under¬ 
taken  to  investigate  this  phenomenon  in  more  detail.  Several  large  events 
from  various  regions  were  selected,  and  signal  and  noise  RMS  levels  were 
computed  for  each  event  on  all  subarray  beams  filtered  with  the  standard 
filter.  The  events  were  grouped  by  region,  and  signal -to -noise  ratios  (SNR) 

(in  dB)  were  averaged  over  all  events  within  each  region  for  each  subarray. 

Since  the  noise  level  does  not  vary  significantly  between  subarrays,  these 
SNR  values  are  essentially  equivalent  to  the  corresponding  signal  ampli¬ 
tudes. 

The  results  are  presented  in  Figures  II - 5  to  II - 8  for  four 
selected  regions:  Kuriles,  Kazakh,  Y unan  and  Kirgiz.  A  map  of  the 
NORSAR  array  is  shown  on  each  figure  together  with  contours  indicating 
the  pattern  of  signal  amplitude  variations  for  each  region.  The  interval 
between  contours  is  2  dB,  and  the  numbers  represent  SNR  level  in  dB 
relative  to  the  average  subarray.  Within  each  region,  typical  standard  de¬ 
viations  for  a  subarray  across  the  ensemble  of  events  were  2  dB,  thus  the 
contour  maps  should  give  a  reasonably  reliable  picture  of  the  NORSAR 
amplitude  distribution  for  the  regions  selected. 
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Region:  Kazakh;  Azimuth:  74.5  -  75.5  degrees;  Distance:  37.7  -  38.  2  degrees 
Events:  KAZ/081/04N  KAZ/115/03N,  KAZ/145/04N,  KAZ/157/04N,  KAZ/181/0'N 

FIGURE  II -8 

CONTOUR  PLOT  SHOWING  SUBARRAY  SNR  VARIATIONS  FOR 
EVENTS  FROM  KAZAKH  REGION 


It  appears  from  these  data  that  the  eastern  part  of  the  NORSAR 
array  (subarrays  01C  through  07C  and  01B  through  04B)  generally  show  higher 
signal  amplitudes  than  subarrays  located  further  west.  This  seems  to  hold 
true,  with  few  exceptions,  for  all  of  Eurasia;  however,  our  limited  data  from 
the  Western  hemisphere  indicates  that  this  pattern  is  not  globally  valid. 

The  great  variation  in  amplitude  patterns  even  for  regions 
fairly  close  together  (like  Kazakh  and  Kirgiz)  is  striking.  Of  particular 
interest  is  the  Kazakh  region,  (Figure  II - 8)  where  subarray  12C  exhibits 
14  dB  lower  SNR  than  the  average  subarray.  This  spread  is  much  greater 
than  has  been  observed  for  other  regions,  including  other  underground  nuclear 
test  sites. 

An  attempt  was  made  to  relate  the  amplitude  variations  to 
the  structure  of  the  Mohorovicic  discontinuity  underneath  the  NORSAR  array. 
Similar  studies  to  explain  amplitude  anomalies  across  the  LASA  array  have  been 
performed  by  Aki  (1973)  and  others. 

Depth  contours  for  the  Moho  underneath  NORSAR  are  presented 
in  Figure  II - 9 ,  as  measured  by  Kanestrom  and  Haugland  (1971).  The  arrow  on 
this  figure  represents  the  wavefront  arrival  for  Eastern  Kazakh  events  (azimuth 
75  degrees),  and  is  seen  to  follow  closely  the  trend  of  a  ridge  in  the  Moho  below 
subarray  05C.  The  arrow  is  furthermore  perpendicular  to  the  depth  contours 
corresponding  to  the  sudden  depth  increase  of  the  Moho  in  western  direction. 

Figures  II - 1 0  and  II -  1 1  are  conceptual  models  of  the  structure 
underneath  the  NORSAR  array  as  seen  by  applying  two  vertical  sections  as  in¬ 
dicated  in  Figure  II - 9 .  The  directions  of  wavefront  approach  indicated  in 
Figures  11-10  and  II - 1 1  correspond  to  events  from  Eastern  Kazakh.  The 
refraction  effects  shown  in  these  figures  seem  to  explain  well,  in  principle, 
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FIGURE  II- 1 1 


CONCEPTUAL  MODEL  (SECTION  LINE  II  IN  FIGURE  II - 9) 
SHOWING  THE  EFFECTS  OF  A  SUDDEN  DEPTH  INCREASE 
OF  THE  MOHO.  SEISMOMETERS  LOCATED  AT  THE 
SHADED  AREA  WILL  SHOW  ABNORMALLY 
LOW  SIGNAL  AMPLITUDES. 
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how  the  seemingly  anomalous  behavior  of  subarray  05C  and  12C  amplitudes 
for  Kazakh  events  may  occur. 

Note  that  Figure  II - 1 0  will  remain  unchanged  for  all  wave- 
fronts  with  the  same  azimuth  as  the  Kazakh  events,  while  the  location  of  the 
shaded  area  in  Figure  II- 11  will  vary  with  the  angle  of  incidence  of  the  wave- 
front.  This  explains  why  Kazakh  is  the  only  region  for  which  extremely  low 
amplitudes  for  subarray  12C  are  seen,  while  low  amplitudes  for  subarray  05C 
are  observed  also  for  events  from  Yunan,  which  also  is  at  75  degree  azimuth 
from  NORSAR. 

It  is  conceivable  that  most  of  the  amplitude  variations  between 
NORSAR  subarrays  may  be  explair  ed  as  caused  by  curvature  of  the  Moho 
boundary  underneath  the  array  together  with  random  effects  due  to  scattering. 
However,  a  much  more  detailed  analysis  will  be  required  to  reach  a  definite 
conclusion  with  respect  to  this  matter. 

An  interesting  consequence  of  the  regional  consistency  in  the 
subarray  amplitude  variations  is  that  diversity  stack  beamforming  with  pre¬ 
defined  weights  would  be  a  realistic  alternative  to  conventional  bcamforming 
for  NORSAR  on-line  detection  processing.  Also,  the  large  spread  in  subarray 
amplitudes  suggests  that  forming  partial  array  beams,  using  only  the  best 
subarrays  for  each  region,  might  improve  signal -to -noi se  ratios  while  reducing 
computational  load  for  the  NORSAR  DP.  This  possibility  will  be  investigated  in 
future  work. 

D.  TIME  DELAY  ANOMALIES 

As  was  stated  in  Special  Report  No.  6,  significant  time  delay 
anomalies  (deviations  from  plane  wave  propagation  along  the  great  circle 
azimuth)  are  observed  for  the  NORSAR  array,  and  must  be  taken  into  account 
in  array  beamforming. 
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Inter -subarray  time  delay  anomalies  were  calculated  for  all 


large  SNR  events  by  computing  the  cro  s  s -co  rr  elation  functions  between  the 
reference  subarray  beam  and  the  remaining  subarray  beams.  A  signal  gate  of 
3  seconds  was  used  in  most  cases.  All  computed  delays  were  hand  checked 
for  inter-subarray  consistency  and  also  compared  with  delays  from  nearby 
events.  Whenever  a  lack  of  consistency  was  found,  the  necessary  adjustments 
were  carried  out  by  the  analyst. 

For  the  teleseismic  regions  (distances  30-100  degrees)  con¬ 
sistent  time  delays  were  generally  obtained  without  difficulty,  thus  confirming 
observations  reported  in  Special  Report  No.  6.  The  problem  mentioned  in  that 
report  with  respect  lo  finding  consistent  time  delays  for  subarray  7  (06B)  for 
Kirgiz  events  is  now  believed  to  have  been  caused  by  the  phase  reversal  in 
sensor  06B03  discussed  in  Section  I  of  this  report. 

As  can  be  expected,  the  irregular  shape  of  the  Mohorovicic 
discontinuity  underneath  the  NORSAR  array  described  in  the  preceding  sub¬ 
section,  has  ramifications  with  respect  to  the  travel  time  anomalies  for  signals 
received  at  NORSAR.  As  an  example,  the  pattern  of  measured  deviations  from 
a  plane  wavefront  for  events  from  Kirgiz  is  shown  in  Figure  11-12.  The  earliest 
signal  arrivals  relative  to  the  plane  wavefront  occur  for  subarrays  01A  and 
08C,  while  subarray  03C  has  a  signal  ai  rival  0.7  seconds  too  late  (relative 
to  the  01A  arrival).  Time  anomalies  of  up  to  half  this  magnitude  may  be 
accounted  for  by  the  relatively  longer  path  travelled  through  the  earth's  crust 
for  rays  arriving  at  03C  than  for  01A.  The  resemblance  of  this  time  delay 
anomaly  pattern  to  the  Moho  pattern  shown  in  Figure  II - 9  indicates  that  the 
remaining  part  of  the  time  delay  anomalies  may  come  from  velocity  structures 
in  the  earth's  mantle  related  to  the  shape  of  the  Moho. 
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FIGURE  11-12 


CONTOUR  PLOT  SHOWING  THE  MEASURED  TIME  DELAY 
DEVIATIONS  FROM  A  PLANE  WAVEFRONT  FOR  EVENTS 
FROM  KIRGIZ.  THE  NUMBERS  REFER  TO  DELAY  IN 
DECISECONDS,  THE  EARLIEST  ARRIVAL  OCCURRING 

AT  0 1 A  AND  08C. 


Near  regional  events,  notably  from  Western  Russia  and  the 
Mediterranean  area,  still  represent  a  considerable  problem  with  respect  to 
finding  appropriate  inter- subarray  time  delays.  However,  the  swarm  of 
events  from  Italy  during  February  4-8,  1972  has  given  some  insight  to  the 
behavior  of  delay  anomalies  for  Mediterranean  region  events. 

For  all  the  Italy  events  processed,  it  was  observed  on  tin- 
time  domain  traces  that  at  least  one  cycle  of  0.  5  Ha  energy  precedes  the  arrival 
of  a  1  Hz  P-wave  (e.  g.  ,  event  ITA/036/01N,  shown  in  Figures  II - 1  1  and  11-14). 
However,  the  relative  sizes  of  the  corresponding  spectral  peaks  were  seen  -Q 
vary  markedly  from  event  to  event.  For  one  event,  ITA/035/02N,  the  0.  5  Hz 
energy  completely  dominates  the  spectrum,  as  seen  in  Figures  11-15  and  11-16. 


For  all  Italy  events  except  ITA/035/02N,  similar  sets  of  time 
delay  anomalies  were  found  applying  the  cro  ss -correlation  procedure.  In¬ 
terestingly,  when  the  standard  filter  was  applied  prior  to  cross-correlat  on 
for  ITA/035/02N  the  resulting  time  delays  were  consistent  with  the  othei  ev'uts. 
(Table  II-l).  Figure  11-17  shows  the  power  spectrum  of  ITA/035/02N  when 
applying  the  delays  computed  from  filtered  data. 

Our  conclusions  from  studying  Italy  events  can  thus  be  stated 

as  follows: 

•  Energy  of  distinctly  different  frequencies  arrives  at  the  NORSAR 

array  at  distinctly  different  times.  This  suggests  the  existence 
of  two  different  propagation  paths  with  highly  different  attenuation 
characteristics  for  Italy  events. 

.  Time  delay  anomalies  computed  for  Italy  events  vary  consider¬ 

ably  with  frequency.  Thus,  array  beamforming  based  on 
delays  computed  from  low  frequency  energy  may  cause  sub¬ 
stantial  loss  for  higher  frequencies 
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FIGURE  11-17 

SUBARRAY  AND  ADJUSTED-DELAY  ARRAY  BEAM  SPECTRA  FOR 
EVENT  ITA/035/02N.  ADJUSTED  DELAYS  IN  THIS  CASE  WERE 
COMPUTED  FROM  FILTERED  SUBARRAY  BEAMS. 
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A  similar  phenomenon  of  different  arrival  times  for  low  and 
high  frequency  energy,  although  much  less  pronounced,  has  been  observed 
for  some  other  events  from  the  Mediterranean  area.  Figure  11-18  shows  as  an 
example  the  waveforms  for  an  earthquake  from  Greece,  GRE/109/02N. 


E.  NORSAR  m^  MEASUREMENTS 

NORSAR  values  were  measured  on  a  number  of  detected 
events  from  Eurasia  using  the  formula: 

m.  =  log  A/  T  +  B 

D 

where:  A  is  the  maximum  peak-to-peak  signal  amplitude  in  m // 

on  the  adjusted  delay  array  beam  (corrected  for  seismometer 
response) . 

T  is  the  period  of  the  cycle  with  the  maximum  amplitude. 

B  is  the  distance  factor. 


Values  for  B  were  given  in  Special  Report  No.  6,  Table  III -4; 

the  same  values  are  also  usea  for  the  calculation  of  LASA  m  's. 

b 

Figures  11-19  and  11-20  are  plots  of  NORSAR  m  's  versus 

b 

either  PDE  m^'s  (dots)  or  LASA  m^'s  (crosses).  Magnitudes  reported  from 
ISM  are  marked  as  circles.  Figure  11-19  comprises  67  events  from  the  Japan 
Kuriles -Kamchatka  region  while  Figure  11-20  is  based  on  142  events  from  the 
remainder  of  Eurasia. 


For  the  Japan  to  Kamchatka  region  NORSAR  m^  values  appear 
to  be  generally  slightly  lower  than  LASA/PDE  values.  The  average  difference 
is  0.  14  magnitude  units,  and  appears  to  be  independent  of  event  magnitude. 
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NORSAR  MAGNITUDE  VERSUS  PDE  OR  LASA  MAGNITUDE, 
KAMCHATKA-KURILES-JAPAN  REGION 


Note  that  most  of  the  low  magnitude  events  in  this  region  have  magnitudes 

from  LASA,  which  is  located  at  approximately  the  same  epicentral  distance 

as  NORSAR.  The  spread  in  the  data  appears  to  be  considerable  in  view  of 

this  fact.  Also,  because  of  the  large  number  of  small  events  not  detected  by 

NORSAR,  a  better  estimate  of  the  bias  is  obtained  by  only  considering 

events  of  PDE/LASA  m  greater  than  4.0;  this  yields  an  average  m  difference 

D  b 

of  0.  20  m,  units, 
b 

All  events  from  other  regions  of  Eurasia  have  been  included 
in  Figure  11-20.  Again,  NORSAR  exhibits  a  generally  lower  m^  (an  average 
difference  of  0.29  magnitude  units)  but  this  difference  seems  to  increase  as 
the  event  magnitude  decreases.  Since  this  did  not  appear  to  be  the  case  for 
the  Japan  to  Kamchatka  region,  it  might  seem  reasonable  to  relate  the  problem 
to  the  following  two  factors: 

•  Several  of  the  low  magnitude  events  in  the  last  case  have  an 

m  from  PDE.  Stations  in  the  PDE  net  that  report  m  for  a 
D  b 

small  earthquake  are  likely  to  have  a  favorable  radiation 
pattern  for  this  event,  and  will  consequently  report  a  re¬ 
latively  high  amplitude. 

•  While  most  of  the  large  events  in  our  population  are  from  a 
distance  of  40  degrees  or  more,  where  beamforming  loss  is 
low,  a  greater  number  of  the  smaller  earthquakes  are  from  the 
Mediterranean  region,  where  large  beamforming  loss  accounts 
for  low  NORSAR  magnitudes. 

However  the  above  problem  still  requires  further  investigation 
before  a  definite  conclusion  can  be  drawn 

Figure  11-21  gives  histograms  of  the  magnitude  difference 
(NORSAR  minus  PDE  or  LASA)  for  the  two  subregions.  The  negative  bias 
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FIGURE  11-21 

HISTOGRAMS  OF  THE  DIFFERENCES  BETWEEN 
NORSAR  AND  PDE/LASA  m  FOR  209  EVENTS 
FROM  EURA&A 


could  be  accounted  for  by  signal  loss  in  array  beamforming,  which  typically 
is  3  to  4  dB  for  teleseismic  events  (0.  15  to  0.  2  magnitude  units).  For  near 
regional  events  the  loss  will  generally  be  higher,  thus  explaining  the  greater 
negative  bias  for  the  continental  Eurasian  region. 

The  distribution  about  the  bias  appears  to  be  approximately 
normal,  and  would  reflect  normal  variation  in  signal  amplitude  due  to  radiation 
patterns  and/or  propagation  effects.  Standard  deviations  in  the  two  event  popu¬ 
lations  are  0.  32  for  the  Japan  to  Kamchatka  region  and  0.  36  for  the  Continental 
Eurasia. 

All  the  NORSAR  magnitudes  utilized  in  the  preceding  analyses 

have  been  measured  by  TI.  For  part  of  the  data,  in  values  were  available 

b 

from  the  NORSAR  seismic  bulletin  compiled  at  Kjeller,  Norway.  A  comparison 

was  carried  out  to  see  if  any  systematic  differences  were  present.  Figure  11-22 

is  a  histogram  of  measured  deviations  for  62  events  from  distances  beyond  30 

degrees.  It  is  seen  that  TI  measurements  have  a  slight  negative  bias  (-0.02  m 

units),  but  that  most  measurements  are  consistent.  The  two  events  that  show 

particularly  large  deviations,  SIB/014/03N  and  L01/058/10N,  have  significant 

location  differences  between  the  NORSAR  bulletin  and  the  PDE  epicenter. 

In  general,  the  observed  differences  are  due  to  different  location  estimates, 

inaccuracies  in  time  delay  corrections  and  lack  of  precision  in  rn  measurements 

b 

It  is  interesting  to  compare  Figure  11-22  to  the  corresponding 
histograms  for  NORSAR  -  PDE/ L ASA  magnitudes  (Figure  11-21).  Figure  11-22 
may  be  viewed  as  expressing  the  typical  uncertainty  involved  in  two  independent 
measurements  of  the  same  mb>  and  the  corresponding  standard  deviation  is 
approximately  0.  15.  This  contrasts  with  a  standard  deviation  of  slightly  above 
0.  30  for  the  interstation  differences  from  Figure  11-21,  thus  confirming  that 
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the  observed  spread  between  NORSAR  and  PDE/LASA  magnitudes  is  mostly  real 
and  not  due  to  random  measuring  effects.  (The  "real"  spread  in  this  case  would 
be  approximately:  \J 0.  342  -  0.  152  =  0.  31.) 

One  systematic  regional  dfference  was  seen  in  the  observed 
data,  the  Iran  region  where  8  events  showed  an  average  negative  bias  of  0.  2 
for  the  TI  measurement.  In  fact  these  8  events  accounted  for  all  of  the  bias 
in  the  total  histogram.  It  thus  appears  that  the  time  delay  corrections  used  by 
TI  were  not  optimum  for  this  region. 

Our  event  population  was  not  sufficiently  large  to  enable  us  to 
carry  out  a  comparison  between  m^  values  measured  by  TI  and  NTNF/NORSAR 
for  near  regional  events,  except  for  the  Italy  region.  For  the  10  events  from 
Italy,  the  TI  measurements  averaged  0.4  mb  units  higher  than  NORSAR's  and 
it  appeared  that  the  NORSAR  Event  Processor  was  operating  with  poor  time 
delay  corrections  for  this  region  at  the  time  these  measurements  took  place. 
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SECTION  III 

ARRAY  PROCESSING  RESULTS 


A.  INTRODUCTION 

This  section  discusses  the  improvement  of  signal -to -noise 
ratio  (SNR)  and  the  signal  degradation  resulting  from: 

•  Subarray  Beamforming 

•  Array  Beamforming 

•  Application  of  the  "standard"  filter  (Figure  III  - 1 ) 

A  total  of  71  events  from  Eurasia  were  used  in  these  analyses; 

the  magnitude  range  of  these  events  was  from  m  =  4.  0  to  6.  0. 

b 

B.  SUBARRAY  BEAMFORMING  PERFORMANCE 

Measurements  of  signal  degradation  and  SNR  improvement 
from  single  sensor  to  subarray  beam  were  made  for  three  events,  WRS/295/05N, 
KAZ/181/04N  and  GRE/074/15N.  The  first  two  of  these  are  presumed  explosions. 
All  three  contain  an  appreciably  greater  than  average  proportion  of  high-frequency 
energy.  Previous  analysis  (Special  Report  No.  6,  1972)  has  shown  that  sub¬ 
array  beams  for  other  events  will  in  general  perform  at  least  as  well  as  those 
discussed  below. 

Signal  degradation  (D)  and  SNR  improvement  (I)  were  computed 
for  each  subarray  using  the  formulas: 

D  =  10  log1()  (S-N)  -  10  log1()  (Sb  -  Nb) 
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FILTER  RESPONSE 


FIGURE  III- 1 

FREQUENCY  RESPONSE  CURVE  OF  THE  STANDARD  FILTER 
USED  FOR  P-WAVE  DETECTION  FOR  NORSAR  DATA 
CORNER  FREQUENCIES  ARE  1. 2  AND  2.  8  Hz 


I  =  10  log 

b10 


10  log 

s10 


N 


where  S  and  Nare  signal  and  noise  powers  over  6.  4  and  70  seconds  windows 

respectively,  averaged  over  all  sensors  within  a  subarray;  while  S  and  N 

b  b 

are  signal  and  noise  powers  for  the  subarray  beam  over  the  same  windows. 

Table  III  - 1  summarizes  subarray  beamforming  results  for 
the  three  events  investigated,  using  both  plane-wave  and  adjusted  delays  for 
subarray  beamforming.  The  adjusted-delay  subarray  beam  shows  an  average 
improvement  of  less  than  0.  5  dB  over  the  plane  wave  beam.  The  gain  result¬ 
ing  from  using  adjusted  delays  in  subarray  beamforming  is  thus  only  marginal 
even  for  close-in,  high  frequency  events. 

An  interesting  observation  is  that  the  average  subarray  beam 
wide  band  noise  reduction,  which  is  the  sum  of  the  SNR  improvement  and  the 
signal  degradation,  is  10.  5,  9  and  11  dB  for  the  three  events  respectively. 

This  is  well  above  the  expected  7.  8  dB  corresponding  to  \/n  noise  reduction. 

The  apparent  explanation  for  this  phenomenon  is  the  strong,  directional  3-6 
second  Rayleigh  wave  noise  field  generally  propagating  from  the  Northwest 
across  the  NORSAR  array.  (Special  Report  No.  5,  1972).  The  typical  sub¬ 
array  geometry  seems  to  take  advantage  of  the  coherency  of  this  noise  for 
beamforming  toward  the  East  and  Southeast.  In  fact,  when  steering  the  sub¬ 
array  beams  toward  the  Northwest  (300  degrees  azimuth,  11  km/sec  velocity) 
the  wide-band  noise  suppression  for  the  noise  preceding  WRS/295/05N  averaged 
only  5.  3  dB  across  the  subarrays. 

C.  ARRAY  BEAMFORMING  PERFORMANCE 

Array  beamforming  SNR  improvement  and  signal  degradation 
were  computed  for  71  Eurasian  events  in  a  manner  similar  to  the  computations 
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SUMMARY  OF  SUBARRAY  BEAMFORMING  ANALYSIS 
RESULTS  FOR  WIDE- BAND  SIGNALS 


performed  in  Subsection  III-B.  Figure  III-2  is  a  histogram  showing  the  dis¬ 
tribution  of  wide -band  SNR  improvement  over  the  entire  event  population.  More 
than  80%  of  the  events  fall  in  the  9-12  dB  bracket;  of  the  six  events  with  gains  less 
than  7  dB,  five  are  from  the  Mediterranean  region  and  one  is  a  presumed  explo¬ 
sion  from  E.  Kazakh  (KAZ/282/06N) . 

Average  SNR  gain  was  10.  1  dB,  while  signal  degradation  over 
the  event  ensemble  averaged  3.  0  dB.  The  sum  of  these  two  numbers  conforms 
well  to  the  expected^N  noise  reduction  across  the  array,  which  amounts  to 
13.  4  dB  for  N  =  22. 

The  array  beamfo rming  gains  were  also  computed  for  each  event 
after  the  standard  filter  had  been  applied  to  the  subarray  beams.  With  this  pro¬ 
cedure,  an  average  SNR  gain  of  9.  3  dB  was  achieved,  while  signal  degradation 
averaged  4.0  dB.  This  higher  degradation  for  filtered  signals  is  to  be  ex¬ 
pected  since  the  filter  attenuates  the  low  frequency  energy,  which  is  generally 
the  most  coherent  part  of  the  signal  across  the  array.  Figure  III  -  3  shows  the 
average  subarray  and  array  beam  spectra  for  a  Turkey  event,  TUR/276/17N, 
and  can  serve  as  a  good  illustration  of  the  abo vementioned  point.  This  event 
has  a  low  array  beamforming  loss  around  1  Hz,  while  beamforming  loss  around 
1.  5  Hz  is  close  to  the  13  dB  noise  suppression. 

Formation  of  dive rsity- stack  array  beams  yielded  an  average 
improvement  over  the  adjusted-dels»y  beam  of  1.0  dB  for  wideband  signals  and 
1. 6  dB  for  signals  to  which  the  standard  filter  had  been  applied.  Table  III-2 
shows  the  distribution  of  the  improvements  across  the  event  ensemble. 

D.  STANDARD  FILTER  PERFORMANCE 

The  SNR  improvement  and  signal  degradation  obtained  by 
applying  the  standard  filter  were  computed  for  the  ensemble  of  71  Eurasian 
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FIGURE  III -2 


GAIN  IN  SIGNAL -TO -NOISE  RATIO  FROM  WIDE -BAND 
SUBARRAY  BEAM  TO  ADJUSTED-DELAY  ARRAY 
BEAM  FOR  71  EURASIAN  EVENTS 


events  previously  referred  to  in  this  section.  Figure  III -4  is  a  histogram 
showing  the  average  subarray  SNR  filtering  gain  for  the  individual  events;  a 
similar  representation  of  SNR  improvements  for  the  filtered  array  beam  is 
presented  in  Figure  III-5. 

Average  SNR  improvement  achieved  by  the  standard  filter 
is  7.  5  dB  on  the  subarray  beam  and  6.  8  dB  on  the  array  beam  level,  thus 
indicating  that  beamforming  tends  to  decrease  the  effectiveness  of  the  standard 
filter.  The  spread  in  this  data  is  considerable,  with  observed  SNR  gains 
ranging  from  -4  to  16  dB  for  the  filtered  array  beam.  In  view  of  the  stable 
spectral  characteristics  of  the  short  period  noise  at  NORSAR  (Special  Report  No.  6 
No.  6,  1972),  these  fluctuations  essentially  reflect  variation  in  signal  spectral 
contents.  E.  g.  for  the  Turkey  region,  with  dominantly  low  frequency  signals 
(Subsection  II-B)  standard  filter  SNR  improvement  averages  only  2  dB,  while 
the  average  gain  for  Afghanistan  events  is  10  dB  and  for  E.  Kazakh  presumed 
explosions  14  dB. 

Signal  power  attenuation  caused  by  the  standard  filter  is  shown 
as  a  histogram  in  Figure  III-6.  Average  signal  loss  is  6.  5  dB  on  the  filtered 
array  beam,  with  values  for  individual  events  ranging  from  1  to  17  dB.  This 
average  loss  combined  with  the  average  SNR  gain  of  6.  8  dB  mentioned  earlier 
accounts  for  an  average  noise  reduction  of  about  13  dB  achieved  by  the  standard 
filter.  This  number  is  consistent  with  the  results  from  our  noise  analysis  pre¬ 
sented  in  Special  Report  No.  6,  1972. 

Figure  III -7  is  a  plot  of  the  effectiveness  of  the  standard  filter 
as  a  function  of  event  wide-band  SNR.  It  appears  that  filter  performance  is 
largely  independent  of  event  size.  There  is  a  slight  trend  toward  higher  filtering 
SNR  gains  for  smaller  events  when  Turkey  events  are  eliminated;  this  might  be 
attributed  to  a  general  tendency  for  smaller  events  to  contain  proportionally 
more  high  frequency  energy  than  large  events.  An  enlarged  data  base  will 
be  necessary  to  study  this  problem  in  more  detail. 
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Signal-to -noise  Ratio  Improvement  (dB) 


FIGURE  III -4 

AVERAGE  SUBARRAY  BEAM  SNR  IMPROVEMENT  ACHIEVED 
BY  APPLICATION  OF  THE  STANDARD  FILTER 
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Average  Subarray  Wide -Band  SNR  (dB) 


FIGURE  III -7 


SNR  IMPROVEMENT  ACHIEVFD  BY  THE  STANDARD  FILTER  ON 
THE  SUBARRAY  LEVEL  AS  A  FUNCTION  OF  EVENT  SNR 
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SECTION  IV 

NORSAR  TELESEISMIC  DETECTION  CAPABILITY 

A.  ESTIMATE  OF  THE  NORSAR  m,  DETECTION  THRESHOLD 

b 

A  study  was  carried  out  to  deter  line  the  detectability  of  tele- 
seismic  P-  waves  using  the  NORSAR  short-period  array.  The  adjusted-delay 
array  beam  was  used  for  P-  wave  detection.  In  almost  all  cases  where  event 
magnitude  was  below  4.  5,  a  "standard"  filter  was  us»d  (See  Figure  III-l). 

The  procedure  leading  to  the  selection  of  this  filter  was  outlined  in  Special 
Report  No.  6 . 

The  data  base  used  in  this  study  consisted  of  a  total  of  344  events. 
Thr  106  events  listed  in  Special  Report  No.  6  were  supplemented  with  238  addi¬ 
tional  events  listed  in  Table  1-1. 

In  order  to  make  possible  an  unbiased  estimation  of  the  NORSAR 
detection  threshold,  all  events  in  this  population  that  had  been  selected  using 
the  NORSAR  seismic  bulletin  as  a  source  were  eliminated.  Also  deleted  were 
events  reported  in  the  bulletin  from  the  International  Seismic  Month  (ISM)  where 
NORSAR  had  been  listed  as  the  primary  source. 

We  thus  arrived  at  a  data  base  consisting  of  303  events  from 
Eurasia.  This  population  should  be  a  representative  set  of  seismic  events  from 
the  Eurasian  Continent,  including  Europe  and  the  Mediterranean  all  of  Continental 
Asia  and  the  Kamchatka-Kuriles -Japan  arc.  The  Philippines  Islands'  region  is  not 
included  and  only  a  few  Japan  events  are  present.  Our  event  population  includes 
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relatively  too  many  high  magnitude  events,  but  this  bias  will  of  course  not 
influence  our  estimate  of  the  NORSAR  incremental  detection  threshold. 

The  data  base  was  divided  into  two  subsets.  121  events  from 
the  Japan-Kuriles-Kamchatka  arc  (called  Japan  region)  and  182  events  from  all 
other  regions  in  Eurasia  (called  Continental  Eurasia).  A  further  regional  sub¬ 
division  would  be  desirable,  but  our  limited  event  population  so  far  does  not 
permit  this. 

Of  the  121  events  in  the  Japan  region  a  total  of  26  events  ranging 
in  body-wave  magnitude  from  3.  3  to  4.  4  were  not  detected  by  NORSAR.  Figure 
IV- 1  is  a  histogram  repr eventing  the  magnitude  distribution  of  Japan  region 
events  processed  and  indicating  how  many  of  these  were  not  detected. 

Of  the  182  events  in  Continental  Eurasia,  a  total  of  12  events 
between  m^  =3.3  and  4.  8  were  not  detected  by  NORSAR.  The  one  event  of 
m^  =4.8  (ITA  035/04N)  not  detected  may  possibly  have  been  assigned  too 
high  a  magnitude,  since  only  one  USCGS  station  reported  a  signal  amplitude 
for  this  event.  This  station  (LOR)  averaged  0.  7  mfa  units  higher  than  the 
PDE  m^  for  other  Italy  events.  Figure  IV-2  is  a  histogram  representing 
the  magnitude  distribution  and  detection  status  for  events  within  Continental 
Eurasia. 

Figure  T.V-3  shows  a  plot  of  the  detected/not  detected  status  for 
events  from  Continental  Eurasia  of  magnitude  less  than  4.  7.  There  is  no  def¬ 
inite  trend  towards  lower  detectability  for  greater  epicentral  distance,  but  our 
population  of  small  events  is  scarce  above  40  degrees  distance.  It  is  still 
interesting  to  note  the  good  detection  performance  between  35  and  55  degrees, 
most  of  the  events  in  this  group  are  from  Central  Asia. 

An  estimate  of  the  incremental  detection  probability  for  the 
NORSAR  array  wa?  derived  from  Figure  IV-4  .  To  stabilize  the  pictures  for 
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FIGURE  IV-3 

NORSAR  P-WAVE  DETECTIONS/MISSES 
VS.  m  AND  EPICENTRAL  DISTANCE 
(EURASIA  APART  FROM  JAPAN -KAMCHATKA) 
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the  two  subregions  individually,  the  events  were  grouped  in  steps  of  0,  2  magm- 
tude  units.  The  detection  probability  in  each  bin  is  the  ratio  of  the  number  of 
events  detected  to  the  total  number  of  events  processed.  It  appears  that  the 
90%  incremental  detection  level  is  around  mb  =  4.  3  for  the  Japan -Kamchatka 
region  and  between  mb  =  4.  1  and  4.  2  for  the  remainder  of  Eurasia.  The  histo- 
gram  representing  all  events  is  somewhat  more  stable,  and  suggests  a  90% 
incremental  detection  threshold  near  =  4.  2  for  all  of  Eurasia. 

The  apparent  discrepancy  between  this  result  and  the  m 
threshold  of  4.  3  -  4.  4  suggested  in  Special  Report  No.  6  seem  to  be  attributed 
to  lack  of  data  and  also  the  fact  that  most  of  the  low  magnitude  events  were 

from  the  Japan-Kuriles  -Kamchatka  region  where  detection  capability  is  some¬ 
what  poorer. 

Our  90%  mb  threshold  of  4.  2  agrees  well  with  the  correspondin' 
LASA  threshold  of  3.9  considering  that  the  NORSAR  RMS  noise  levels  through 
the  standard  filter  is  about  a  factor  of  2  higher  than  those  observed  at  LASA 
(Special  Report  No.  6);  (Dean  1971).  This  accounts  for  a  difference  of  about 
0.  3  mb  units  in  detection  threshold  between  the  two  arrays. 

Some  qualifications  still  exist  with  respect  to  the  above  estimate: 
First  of  all,  the  event  population  is  still  scarce,  and  should  be  expanded  further. 
Notably  more  events  of  magnitudes  4.  0  -  4.  5  are  needed.  Secondly,  all  decision, 
regarding  detection/no  detection  have  been  made  by  an  analyst  looking  for  a 
signal  with  a  known  source  location  and  arrival  time,  thus  the  degradation  in¬ 
herent  in  an  automatic  detection  system  has  not  been  a  factor  here. 


B. 


COMPARISON  OF  TI  AND  NORSAR  DETECTION  RATES 


For  the  time  period  I  January  to  20  March  1972  a  study  was 
carried  out  to  compare  the  number  of  events  reported  in  the  NORSAR  seismic 
bulletin  to  the  number  of  detections  claimed  by  the  TI  analyst.  A  total  of  122 
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events  reported  from  PDE/LASA/ISM  were  available  for  this  time  period,  53 
from  the  Japan  to  Kamchatka  region  and  69  from  other  parts  of  Eurasia.  The 
results  are  presented  in  Table  IV- 1. 

It  appears  that  the  NORSAR  bulletin  reported  almost  all  de¬ 
tectable  events  of  magnitude  greater  than  4.  2  for  the  Japan  region  ,  while  several 
events  up  to  magnitude  4.  6  from  other  parts  of  Eurasia  were  missed.  Some  re¬ 
servations  must  be  made  due  to  the  scarcity  of  events  fromthe  Japan  region,  but 
our  data  still  reflects  the  excellent  beam  coverage  of  this  area  for  the  NORSAR 
on-line  Detection  Processor. 

Several  low  magnitude  events  have  not  been  reported  in  the 
NORSAR  bulletin,  this  is  in  part  due  to  a  conservative  acceptance  threshold 
for  the  NORSAR  Event  Processor  at  the  time,  but  may  also  reflect  the  inherent 
limitations  of  an  automatic  event  detector. 

Based  on  our  analysis  of  NORSAR  data  so  far,  it  is  possible 
to  establish  a  simple  theoretical  model  for  the  NORSAR  detection  capability 
and  to  give  tentative  estimates  of  the  parameters  as  follows: 

For  an  arbitrary  earthquake  of  PDE/LASA  body-wave  magni¬ 
tude  m,  NORSAR  can  be  expected  to  see  this  signal  as  a  magnitude  m^  event 

where  m  is  sampled  from  a  Gaussian  distribution  with  mean  m  -  A  and  standard 
N 

deviation  cr  .  (e.  g.  for  the  Japan  to  Kamchatka  region,  0.2  and  o\~  0.3 

according  to  subsection  II-E). 

Whether  or  not  an  event  is  detected  at  NORSAR  is  dependent 
upon  noise  level  and  processing  losses.  In  Special  Report  No.  6  the  average 
RMS  filtered  array  beam  noise  was  found  to  be  0.  12  mfi  ,  and  a  tentative 
standard  deviation  is  2  dB.  Signal  attenuation  by  the  standard  filter  varies  con¬ 
siderably;  values  from  1  to  7  dB  are  typical  for  teleseismic  region.  If  signal- 
to-noise  ratio  (SNR)  required  for  detection  is  set  at  11  dB  (which  is  consistent 
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TABLE  IV-1 


COMPARISON  OF  THE  NORSAR  SEISMIC  BULLETIN  TO  EVENTS 
VERIFIED  BY  TI  USING  NORSAR  DATA 


Magnitude 

Bracket 

Japan -Kamchatk 

a  Arc 

Remainder  of  Eurasia 

Total 

Events 

TI 

Detection 

NORSAR 

Bulletin 

Total 

Events 

TI 

Detection 

NORSAR 

Bulletin 

3.  3  -  3.4 

1 

0 

0 

1 

o 

0 

3.5  -  3.6 

6 

2 

0 

2 

1 

0 

3.  7  -  3.  8 

11 

8 

3 

3 

3 

0 

3.  9  -  4.0 

10 

6 

3 

6 

6 

3 

4.  I  -  4. 2 

8 

6 

4 

9 

8 

3 

4.  3  -  4.4 

3 

3 

3 

8 

•  8 

4 

4.  5  -  4.  6 

5 

5 

4 

12 

11 

7 

4.  7  -  4. 8 

5 

5 

5 

12 

11 

11 

4.9  - 

4 

4 

4 

16 

16 

l 

15 
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with  the  NORSAR  detection  processor),  we  find  that  the  smallest  signal 
generally  detectable  on  an  unfiltered  array  beam  has  a  RMS  value  of  0.  7  mu 
or  a  peak-to-peak  amplitude  of  2  m fi  .  At  a  65  degrees  range  this  corresponds 
to  an  earthquake  of  NORSAR  magnitude  of  approximately  3.75  if  a  period  of 
1.0  seconds  is  assumed.  This  limiting  magnitude  m  varies  with  background 

J-J 

noise  and  as  a  function  of  system  loss,  and  we  will  assume  that  it  follows 
a  Gaussian  distribution  with  expectation  and  standard  deviation  cr 
Values  of  m^~  3.75  for  the  Japan-Kamchatka  region  and  cr  ^  zr  0.  2  are  suggested 
by  the  preceding  considerations. 

We  are  now  able  to  formulate  the  probability  of  NORSAR 
detecting  an  event  of  PDE/LASA  magnitude  m  as  follows: 


m  -  A  -  m 


Pr  (Detect  m)  =  Pr  (m^  >  m^)  =  4>  (  — 


(X 


2  2  2 

where  cr  =  cr  ^  +  cr  ^ 


since  (m  -  m  )  is  Gaussian  with  expectation  (m  -  A  -  m  )  and  standard 
N  L  t) 

deviation  cr  .  <J)  denotes  here  the  cumulative  normal  distribution  function. 

The  detectability  curve  with  the  above  parameters  for  the 
Japan-Kamchatka  area  is  sketched  in  Figure  IV-5.  It  is  seen  to  agree  rea¬ 
sonably  well  with  the  observed  detection  rates  shown  in  Figure  IV-4,  although 
the  90%  detection  level  is  slightly  higher  (around  4.  4)  for  the  theoretical  curve. 

The  above  model  may  be  used  for  estimating  the  total  number 
of  detections  N  by  the  NORSAR  array  for  a  given  seismic  region,  using  the 
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FIGURE  IV- 5 


EXPECTED  DETECTION  CAPABILITY  OF  THE  NORSAR  ARRAY 
FOR  THE  JAPAN-TO-KAMCHATKA  REGION  BASED  ON 
MEASURED  SEISMIC  NOISE  LEVEL 
AND  PROCESSING  LOSSES 


standard  seismicity  formula  for  the  cumulative  number  N  of  earthquakes 
above  a  given  magnitude  m: 

log  N  =  a  -  b  m 
c 

N  is  then  determined  by 


dN  m-A-m 

/oo  p  n 

- —  $  ( - )  dm 

dm  tr 

- 00 

which  yields 

b2  2 

log  *  =  a  -  b  {  mQ  +  A  )  +  - - -  In  10 

Values  cf  a  =  6.  85  and  b  =  1.  04  for  the  year  1965  have  been 
obtained  by  Evernden  (1970)  for  the  Kamchatka/Kuriles  region  of  the  USSR. 
Based  on  this  our  estimate  of  the  NORSAR  optimum  detection  capability  for 
this  area  becomes  approximately  800  events  annually.  Of  course  some  uncer- 
tainty  in  this  number  is  caused  by  variable  seismicity  from  year  to  year. 

Apart  .rom  this,  one  significant  possible  error  source  is  our  estimate  of  the 
50%  incre  nental  detection  threshold  mQ  +  A  .  An  error  of  0.  1  in  this  estimate 
will  cause  our  estimate  of  N  to  be  wrong  by  20-25%. 

The  model  developed  here  does  not  apply  directly  to  the 
general  Eurasian  region  since  distance  factors  vary  greatly  when  near  re¬ 
gional  events  are  included  and  also  the  NORSAR  magnitude  bias  seems  to  show 
a  significant  regional  dependence.  More  data  will  be  needed  to  analyze  theo¬ 
retical  detection  capability  for  the  remainder  of  Eurasia  on  a  regional  basis. 
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SECTION  V 

SHORT  -  PERIOD  DISCRIMINATION 


A.  DEFINITION  OF  DISCRIMINANTS 

1.  P30  Mean  Square 

This  discriminant,  which  is  a  measure  of  event  complexity,  is 
computed  by  crosscorrelating  4  sec  of  the  waveform  (beginning  a  few  points 
before  F-wave  onset)  with  the  next  30  seconds  of  the  waveform  and  with  the 
noise  preceding  the  signal.  A  mean  square,  weighted  by  the  lag,  is  then  com¬ 
puted  from  the  correlations  over  both  30  seconds  of  the  noise  and  30  seconds 
of  the  signal.  The  noise  mean  square  is  subtracted  from  the  signal  mean  square 
to  obtain  the  discriminant  used  (Texas  Instruments  Incorporated,  1971). 

2.  Autocorrelation  Mean  Square 

This  discriminant  is  also  a  measure  of  complexity.  The  auto¬ 
correlations  of  a  30-second  noise  gate  and  of  a  30-second  signal  gate  are  computed 
and  a  weighted  mean  square  then  derived  from  these  correlations  for  the  noise 
and  signal.  The  discriminant  is  derived  from  the  signal  mean  square  minus  the 
noise  mean. 

3.  Envelope  Difference 

This  discriminant  is  also  derived  from  the  P30  correlation  by 
computing  the  mean-square  difference  between  the  envelope  correlation  and  a 
fixed  decaying  exponential,  the  decay  rate  of  which  is  the  average  rate  for  an 
ensemble  of  16  explosions  recorded  at  LASA.  As  with  the  first  two  statistics, 
envelope  difference  is  a  measure  of  complexity. 
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4. 


Dominant  Period 


This  discriminant  is  computed  by  finding  the  cycle  in  the  wave¬ 
form  with  the  maximum  absolute  amplitude;  the  dominant  period  is  the  duration 
of  this  cycle  in  seconds.  This  parameter  can  be  estimated  with  some  confidence, 
even  for  events  with  a  relatively  low  signal-to -noise  ratio,  The  dominant  period 
discriminant  is  a  rough  measure  of  spectral  energy  distribution. 

5.  Spectral  Ratio 

This  discriminant  is  derived  from  the  signal  power  spectrum 
over  a  gate  beginning  just  before  the  signal  arrival.  The  power  spectrum  is 
smoothed  over  three  frequency  points,  and  the  power  in  three  bands  is  com¬ 
puted;  Band  1  :  0.  -  0.  55  Hz;  Band  2  :  0.  55  -  1.  5  Hz;  Band  3:  1.5  -5.0  Hz. 

These  bands  have  been  selected  based  on  NORSAR  data.  Spectral  ratios 
computed  were  Band  3  to  Band  2  and  Band  3  to  Band  1  respectively. 

The  spectral  ratio  that  seems  to  produce  the  best  separation 
for  L ASA  data,  (0.  35-0.  85Hz  to  1 . 45-1.  85Hz)  also  was  computed  for  each  event, 
but  did  not  appear  to  be  as  effective  for  NORSAR  data. 

B.  NORSAR  SHORT-PERIOD  DISCRIMINATION  RESULTS 

Short  period  discriminant  values  for  the  discriminants  defined 
in  V  A  are  plotted  as  a  function  of  body-wave  magnitude  for  a  total  of  269  events 
in  Figures  V-l  to  V-12.  Shallow  earthquakes  and  earthquakes  of  unknown  depth 
are  represented  by  a  cross.  Deep  earthquakes  (of  depths  greater  than  100  km) 
are  denoted  by  a  triangle.  Presumed  explosions  are  indicated  by  an  asterisk. 
Events  from  the  Western  Hemisphere  are  surrounded  by  a  circle. 


Reference  Log  P30  Mean  Square 


2.0  A 


1.5  -j 

J 


1.0  -I 


0,5  4 


0,  ]- 


-0,5  A 


-1.0 


+ 

4- 


+ 

+ 

t  * 


+  +  *  T+ 

+  +++  *  **+  $ 


*++ 

J.  *f 

+A+ 


& 


+ 

+  $ 


+*J++++*£+,  4  + 

*  +  .*  +4.  + 

*  I  +  * 

+  *  ^  a  #  f  + 


+ 

+ 

+  + 
+ 


+ 


4.  J.  zy  , 

1+'^+  *+* 
+  +  +  +  ,  +4. 

+  4-  Af  ±  ++ 


*€? 


+  *-  I*  tT+  + 


+  4-  +  *+  +  4. 

:♦  +  +  +:*v 

+*  ±i*+++j 

+x+A+  *  +  + 


,  ,  +  i'4A+  +  .  . 

++  *  +  *  ~ 

+  $  +  -4® 

,  ,  +A  + 


Hi 


+  A 
+ 


'X 


X  * 


» 


XX 


"1 


,  _  magnitude 

t  Derr  Earlhquake  <shallo»  or  unknown  depth) 
.  ep  Eurastan  Earthquake  (100  km  or  greater) 

*  Presumed  Explosion  from  Eurasia 
North  American  Earthquake 
Presumed  Explosion  from  North  America 
FIGURE  V-  1 


heeerence  subARRAY  P3„  mean  souare  discriminant 


V-3 


. .  - . -  - 


hi  huAuShSBttattM,, 


Square 


ean  Square 


+ 


CL  H 


S 

c 

o 

■  rH 
4-* 

rt 

r-H 

<L> 

U 

o 

u 

o 

4-4 

d 


GO 

o 

-I 

£ 

cti 

<U 


,0 


+i  * 

+  4.  + 


* 


®  t 


4* 


+  +  ++++ 


*V  + 


+ 

*  *t 


++  +¥ 
&  + 


© 


+'*f  :**  tii  :: 

\i*+:  ^  +;+>:+ 


++  4* 


+ + 


+ 

+ 


t 


++i+' , .* 


+  '*’±4.+  +  + 

!*♦+&++  + 

+  A  .•  + 


++++ 

JK 


t :  *v 

+  +4.+  +  +  ±* 

4.  +++a  ++  i 


* 


++ 

A 

+ 


♦x  +' 

+a 


+  3K 


» 


* 

* 


* 


* 


•t-  A  4 

+  A 

+  A 

+  ® 

+  + 


JK 


® 


+  +  A 


',TTi  nTrrrriTrnT 


u  (J 


3.5 


_  ,  ■  ■  T  l-TTI  I  |  |~| 

^  0  '{.  5  5.  0 

MAGNITUDE 


5.5 


nrrrn 


+ 

A 

* 


© 

© 


Eurasian  Earthquake  (shallow  or  unknown  depth) 
Deep  Eurasian  Earthquake  (100  km  or  greater) 
Presumed  Explosion  from  Eurasia 
North  American  Earthquake 
Presumed  Explosion  from  North  America 


FIGURE  V-4 

ADJUSTED-DELAY  BEAM  AUTOCORRELATION  MEAN  SOUARE 
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ADJUSTED-DELAY  BEAM  DOMINANT  PERIOD  (SECONDS) 
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As  was  observed  in  Special  Report  No.  6,  the  discriminants 
appear  to  be  not  very  effective  lor  deep  earthquakes.  In  particular  this  applies 
to  the  autocorrelation  and  envelope  difference  criteria,  while  the  spectral  ratio 
discriminants  seem  to  perform  better  in  this  case.  Also,  it  is  evident  that 
presumed  explosions  from  the  Western  Hemisphere  do  not  follow  the  same 
palteins  as  Eurasian  ones,  thus  we  will  confine  the  remainder  of  this  section 
to  the  discussion  of  NORSAR  short  period  discriminants  applied  to  Eurasian 
events. 


A  comparison  between  the  corresponding  plots  for  the  dis¬ 
criminants  applied  to  the  reference  subarray  beam  versus  the  adjusted  delay 
array  beam  reveals  that  the  array  beam  yields  consistently  as  good  or  better 
separation  between  Eurasian  earthquakes  and  presumed  explosions.  Table  V-l 
was  compiled  by  visually  inspecting  the  plots  shown  and  listing  all  "difficult" 
events  of  magnitude'  5.0  and  greater,  i.  e.,  events  that  would  appear  to  be  within 
the  wrong  population  by  at  least  one  of  the  criteria  applied  to  the  array  beam. 

It  is  seen  from  Table  V-I  that  the  classification  of  an  event 
may  vary  considerably  from  one  discriminant  to  the  other  .  The  presumed 
explosions  that  appeared  to  be  most  difficult  to  classify  were  from  Western 
Russia  WRS/277/10N  and  WRS/295/10N.  These  two  events  Ixrth  had  complex 
waveforms  and  significant  energy  content  below  1.  5  Hz.  Only  the  spectral 
ratio  criter!on  that  compared  high  frequency  energy  to  very  low  frequency 
energy  was  able  to  classify  these  events  properly. 

One  deep  earthquake,  KHR/099/15N,  consistently  behaved  as 
an  explosion  with  respect  to  all  the  discriminants.  One  other  deep  earthquake, 
rAL  /  27 1  /  1  4N ,  and  one  shallow  *a  r  thcpiak  e  KM  1/075/12N,  were  misclassified 
by  the  signal  complexity  c  riteria,  but  behaved  normally  with  respect  to  the 
spectral  content  discriminants.  The  three  remaining  "difficult"  earthquakes 
were  misclassified  by  only  one  of  the  six  discriminants. 
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EFFECTIVENESS  OF  VARIOUS  SP  DISCRIMINATION  CRITERIA 
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It  appears  from  our  da‘a  that  the  discriminants  based  on  spectral 
contents  generally  separate  better  between  earthquakes  and  presumed  explosions 
than  the  discriminants  based  on  waveform  complexity.  However  some  kind  of 
a  multivariate  criterion  would  probably  give  the  most  reliable  performance. 

For  example,  selecting  all  events  that  definitely  belong  to  the  presumed  ex¬ 
plosion  population  for  at  least  one  of  the  spectral  contents  discrim  inants 
(code  E  +  in  Table  V-l)  would  cause  all  the  14  Eurasian  presumed  explosions 
to  be  included,  while  only  1  of  the  51  Eurasian  earthquakes  of  magnitude  5.  0 
and  greater  would  satisfy  the  criterion.  This  one  earthquake  is  KUR/099/15N 
of  depth  126km. 

More  data  from  presumed  explosions  of  magnitudes  below  5.  0 
is  needed  to  determine  the  capabilities  of  short  period  discriminants  applied 
to  low  magnitude  events  and  to  determine  what  kind  of  multivariate  discriminants 
would  be  most  effective  on  NORSAR  data. 


A  comparison  between  the  performance  of  these  Short  Period 

discriminants  and  the  three  combined  Short  Period-Long  Period  discriminants 

(Ms  -  m  ,  AR  -  m  ,  AL  -  m  )  described  in  Special  Report  No.  5  was  carried 
b  b  b 

out.  However,  only  49  events,  8  of  which  arc  presumed  explosions,  have  been 
processed  in  common  for  the  SP  and  LP  evaluation  programs,  thus  our  present 
data  is  sufficient  only  to  yield  a  tentative  indication  cf  the  relative  efficiency 
of  the  two  types  of  criteria.  Of  the  eleven  "difficult"  events  listed  in  Table  3, 
three  had  been  selected  for  LP  discrimination;  KAZ/356/06N,  TIB/123/00N 
and  KUR/213/02N.  All  of  these  events  were  very  clearly  assigned  to  the  proper 
category  by  the  three  SP-LP  discriminants.  In  fact  no  event  in  our  limited 
common  population  failed  any  of  the  combined  SP-LP  discrimination  criteria, 
but  one  event  appeared  to  be  a  borderline  decision,  BLS/210/19N.  This  pre¬ 
sumed  earthquake  of  m  =4.5  had  a  surface  wave  magnitude  of  M  =  2.  5 

b  s 


V-l  7 


However,  all  our  short  period  discriminants  placed  this  event  definitely  in 
the  earthquake  population. 

As  expected  these  preliminary  data  indicate  that  the  Lon« 
Period  discriminants  yield  better  separation  between  oarthquaxes  and  presumed 
explosions  than  our  Short  Period  criteria.  The  questio  of  to  what  degree  SP 
discriminants  can  supplement  LP  criteria  when  those  fail  cannot  be  fully 
answered  from  our  present  data  base.  Future  efforts  will  be  directid  towards 
increasing  the  number  of  common  events  for  and  LP  discrimination  as  well 
as  extending  the  total  data  base  for  discrimination  studies. 
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SECTION  VI 

CONCLUSIONS  AND  FUTURE  PLANS 


A.  CONCLUSIONS 

Conclusions  about  the  performance  of  the  short  period  NORSAR 
array,  based  on  analysis  of  more  than  300  signals  (primarily  from  Eurasia)  are 
given  below. 

Data  quality  is  excellent.  For  about  one-half  of  the  samples 
all  132  sensors  were  operational.  The  worst  data  loss  encountered  was  24 
sensors.  For  six  events  spikes  were  found  in  the  data,  but  these  events  could 
still  be  processed.  Phase  reversals  were  observed  for  a  few  seismometers 
during  parts  of  1971,  but  this  problem  had  been  corrected  for  1972  data. 

Major  conclusions  from  the  signal  analysis  are: 

»  Amplitude  variations  across  the  subarrays  arc  generally  large 

(typically  4:1)  and  show  strong  dependance  upon  source  region. 
However,  within  narrow  regions  a  high  degree  of  consistency 
is  seen,  and  it  appears  that  most  of  the  amplitude  variations 
may  be  explained  by  scattering  effects  due  to  the  irregular 
structure  of  the  Mohorovicic  discontinuity  underneath  the 
NORSAR  array. 

•  lime  delay  anomalies  (deviation  from  plane  wave  propagation 

along  the  great  circle  path)  are  significant  between  subarrays, 
but  a  consistent  set  of  anomalies  can  in  general  be  obtained 
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for  events  of  epicentral  distance  greater  than  30  degrees  from 
NORSAR.  Part  of  the  observed  anomalies,  but  not  all,  can 
be  explained  from  the  variations  in  depth  of  the  Moho  below 
NORSAR. 

For  certain  close-in  events,  notably  events  from  Italy,  low 
and  high  frequency  signal  energy  appear  to  follow  different 
paths;  the  0.  5  Hz  energy  arriving  at  NORSAR  1  -  2  seconds 
before  the  1  Hz  energy.  Significant  differences  were  found 
in  this  case  between  time  delay  anomalies  for  the  low  and  high 
frequency  signal  bands. 

Considerable  variation  in  regional  signal  characteristics  has 
been  observed.  Signal  waveform  complexity  showed  the  expected 
dec  real  e  with  increasing  distance,  except  for  some  high  com¬ 
plexity  events  from  Taiwan  and  South  Kamchatka.  Signal  spectral 
contents  was  more  unpredictable,  with  low  frequency  signals 
being  observed  mainly  from  Italy,  Turkey,  Kirgiz  and  Taiwan 
while  Greece,  Tadzhik  and  Kurile  Island  events  generally 
produced  high  frequency  signals. 

Significant  spread  has  been  observed  between  LASA/PDE  and 
NORSAR  magnitudes  (as  measured  by  the  TI  analyst).  NORSAR 
magnitudes  are  generally  lower,  with  an  average  negative  bias 

of  0.  2  -0.3  mb  units,  and  with  a  standard  deviation  of  0.  3  around 
this  bias. 


The  following  conclusions  werederived 

ing  performance: 


concerning  array  process- 


Average  wide  band  subarray  beam-to-array  beam  SNR  improve 
ment  was  10  dB,  with  80  percent  of  all  examined  events  having 
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values  between  9  .  nd  12  dU.  As  expected,  array  Bain  is  con¬ 
siderably  lower  for  close-in  events  and  presumed  explosions 
with  high  dominant  frequency. 

Average  wide-band  signal  degradation  from  subarr^y  to  array 
level  is  3  dB.  Together  with  the  10  dB  SNR  improvement,  this 
is  consistent  with  the  expected  13.4dB  noise  reduction  for  a  22- 
<?iement  array. 

Diversity  stack  array  beamforming  yielded  an  average  improve¬ 
ment  over  the  adjusted  -  delay  array  beam  of  1.0  dB  for  un¬ 
filtered  and  i  6  dB  for  filtered  signals  (standard  filter). 

Compensating  for  time  delay  anomalies  on  the  subarray  level 
yielded  an  average  of  only  0.  5  dB  SNR  improvement  for  thr-e 
close-in  events  compared  with  plane  wave  subarray  beams. 

Thus  plane  wave  delays  appear  to  be  adequate  on  the  subarray 
level. 

SNR  improvement  achieved  with  the  standard  filter  averaged 
7dB,  with  values  for  individual  events  ranging  from  -4  to  16  dB. 
Significant  regional  dependence  was  seen  m  these  numbers. 

Filter  signal  suppression  averaged  around  6  dB,  and  also  showed 
a  large  variability  between  event?. 

Our  conclusions  concerning  NORSAR  detectability  are: 

90%  incremental  detection  threshold  is  close  to  4.  2  for  all 
of  Eurasia  combined,  and  slightly  nigher  (4.  3)  for  the  Japan  to 
Kamchatka  arc. 

The  number  of  events  reported  in  the  NORSAR  seismic  bulletin 
for  Janua  -y-March  1972  appears  to  be  considerably  lower  than 
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would  be  expected  from  our  analysis  of  NORSAR  detectability. 
Presumably,  this  is  due  partly  to  conservative  prethreshold 
for  operation  of  the  NORSAR  Event  Processor,  partly  to  the 
inherent  limitations  of  the  automatic  signal  detector. 

A  theo  1  etical  model  based  on  'JORSAR  seismic  noise  level 
and  processing  losses  seems  to  give  detectability  estimates 
which  are  cons,  ’tent  with  our  experimental  results. 

Conclusions  with  respect  to  short  period  discrimination  are: 

Discriminants  based  on  spectral  energy  distribution  .' e e m  to 
be  superior  to  discriminan  ^  based  on  the  complexity  af  the 
signal  waveform. 

No  single  discriminant  was  able  to  separate  completely  be¬ 
tween  presumed  explosions  and  earthquakes.  The  best  separa¬ 
tion  was  obtained  by  considering  the  spectral  ratio  of  energy  in 
the  bands  1  5  to  5.0  Hz  and  0  -  0.  55Hx,  although  reservations 
must  be  taken  due  to  possible  bias  caused  by  the  high  signal - 
to-noisc  ra.ios  for  all  events  in  the  presumed  explosion  popu¬ 
lation. 

It  appeared  possible  to  improve  separation  significantly  by 
considering  a  combination  of  discriminants.  However,  no 
formal  multivariate  model  was  established  to  determine  an 
optimum  criterion. 

A  preliminary  study  of  the  performance  of  short  period  dis¬ 
criminants  versus  that  of  M  -  m,  ana  other  SP-LP  discri  - 

s  b 

minants  gave  the  expected  result  that  the  latter  ones  in  general 
produce  a  better  separation  between  earthquakes  and  presumed 
explosions. 
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B.  FUTURE  PLANS 


Future  NORSAR  short  period  evaluation  efforts  will  concentrate 
on  increasing  the  ensemble  of  low  and  intermediate  magnitude  Eurasian  events 
As  the  event  population  increases,  more  emphasis  will  be  laid  upon  obtaining 
regional  estimates  of  signal  processing  losses,  time  delay  and  amplitude  anoma¬ 
lies  and  detection  thresholds.  Regional  signal  characteristics  will  be  investi¬ 
gated  in  furtier  detail. 

The  effects  on  array  gain  from  eliminating  subarrays  that  give 
consistently  low  signal  amplitudes  for  certain  regions  will  he  investigated. 

Measurements  of  short  period  discriminants  will  be  updated, 
and  more  investigation  will  be  performed  to  find  the  optimum  bands  for  com¬ 
puting  spectral  ratio.  The  possibilities  of  multivariate  short  period  discri¬ 
mination  will  he  studied  as  well  as  the  combination  of  short  and  long  period 
discriminants. 
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